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ABSTRACT 



Aims. We aim to investigate the ability of simple spectral models to describe the GRB early afterglow emission. 
Methods. We performed a time resolved spectral analysis of a bright GRB sample detected by the Swift Burst Alert Telescope and 
promptly observed by the Swift X-ray Telescope,with spectroscopically measured redshift in the period April 2005 - January 2007. 
The sample consists of 22 GRBs and a total of 214 spectra. We restricted our analysis to the softest spectra sub-sample which consists 
of 13 spectra with photon index > 3. 

Results. In this sample we found that four spectra, belonging to GRB060502A, GRB060729, GRB060904B, GRB061 1 lOA prompt- 
afterglow transition phase, cannot be modeled neither by a single power law nor by the Band model. Instead we find that the data 
present high energy (> 3 keV, in the observer frame) excesses with respect to these models. We estimated the joint statistical signif- 
icance of these excesses at the level of 4.3 <r. In all four cases, the deviations can be modeled well by adding either a second power 
law or a blackbody component to the usual synchrotron power law spectrum. The additional power law would be explained by the 
emerging of the afterglow, while the blackbody could be interpreted as the photospheric emission from X-ray flares or as the shock 
breakout emission. In one case these models leave a l.lcr excess which can be fit by a Gaussian line at the energy the highly ionized 
Nickel recombination. 

Conclusions. Although the data do not allow an unequivocal interpretation, the importance of this analysis consists in the fact that we 
show that a simple power law model or a Band model are insufficient to describe the X-ray spectra of a small homogeneous sample 
of GRBs at the end of their prompt phase. 

Key words. GRB: X-ray afterglow 



1. Introduction 



, The X-ray telescope (XRT, Burrows et al., 2005l l on board the 
■ Swift satellite jGehrels et al. 2004l i aUows to perform time re- 
' solved spectroscopy of a large number of gamma ray burst 
(GRB) afterglows in the 0.3-10 keV energy band. The vast ma- 
jority of the spectra can be modeled by a single absorbed power 
, law (SPL) model. In most of the afterglows a strong spectral 
1 evolution is observed in the early phases with spectral indexes 
varying in the ra nge 0.5-5 jO'Brien et al. 2006IButler 2007al 
[Zhang et al. 2007 1. These observations are consistent with the 
classical fireball model which describes the burst and afterglow 
emission due to synchrotron radiation. It is possible that detec- 
tion of different spectral features, signature of other emission 
mechanisms, like thermal components or recombination lines 
might give some insight about GRB progenitors, chemical com- 
position, physical conditions and geometry of the GRB environ- 
ment. Before the Swift mission, several X-ray line detections 
were reported in GRB late (> 10 hr) afterglow observations. The 
statistical significance of these detections has been questioned 
by Sako et al. (2005)| who concluded that there were no credi- 
ble X-ray features in any GRB afterglow. Butler et al. (2005), 
however, showed that for GRBOl 1211 the different estimates of 
the statistical significance can be explained by the different ap- 
proach in the continuum modelling. 
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In the Swift afterglow observations only a few deviations 
from SPL have been reported. Most of them have been explained 
by the curvature of the synchrotron spectrum and the presence of 
the vFy peak (Epeak) within the XRT band (Falcone et al. 2006, 
Butler & Kocevski 2007, Goad et al. 2007, Mangano et al. 2007, 
Godet 2007a). 

Moreover Butler (2007) found anomalous soft X-ray emis- 
sion in the spectra of four GRBs which can be interpreted as 
multiple emission lines due to K shell transition in light metals 
as well as thermal emission from a blackbody with temperature 
~0. 1 keV. Gr upe et al. (2 007) and Godet et al. (2007a) found that 
early afterglow data of GRB060729 and GRB050822, respec- 
tively, can be fitted with a SPL plus a blackbody with decreasing 
temperature in the first few hundreds seconds from the beginning 
of the prompt emission. Campana et al. (2006) found a cooling 
thermal component in the spectrum of SN2006aj/XRF060218, 
interpreted as due to the supernova shock break-out (this in- 
terpretation has been subsequently questioned by Li 2007 and 
Ghisellini et al. 2007). 

Starting from the idea that any deviation from a SPL spec- 
tral model, if present, would be a faint signal mostly covered 
by the high level "noise" of synchrotron emission, we searched 
for the best observational conditions to detect it. In particular, 
we searched for high energy excesses with respect to the SPL 
when the spectrum is steepest and, at least in the hard part of the 
energy band, could be the non dominating component. 
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Fig. 1. The 13 very soft spectra are shown together with their best SPL fit. They are sorted in time from the top-left comer. The 
observed time interval and the photon index are indicated in the figures. The asterisks denote the four spectra for which we found 
significant departures from SPL model. 
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Throughout this paper, all errors are quoted at 68% confi- 
dence level for one parameter of interest, unless otherwise spec- 
ified. The reduced will be denoted as x\ ^"d the number of 
degrees of freedom with the abbreviation "d.o.f.". We follow the 
convention Fy{v, t) cc r"v^^, where a and /3 are the temporal 
decay slope and the spectral index, respectively. As time ori- 
gin, we will adopt the Burst Alert Telescope (BAT) trigger (TO). 
The photon index is F = I + /3. Last, we adopt the standard 
"concordance" cosmology parameters, = 0.27, Oa - 0.73, 
ho = 0.71. 

2. Data analysis and sample selection 

We considered the GRB sample detected by BAT, promptly ob- 
served by XRT in the period April 2005-January 2007 which 
collected at least 800 photons. We restricted our analysis to the 
sample with spectroscopically measured redshift in order to sep- 
arate the local contribution to the total absorbing column from 
the Galactic one. We excluded GRB060218 from our analysis 
because of its peculiarity (Campana et al. 2006). For each burst, 
when possible, we split the XRT data in different time intervals 
in such a way that in each of them there are 2000 photons, be- 
fore the background subtraction (the last spectrum collects the 
remaining photons). The final sample consists of 22 bursts for 
a total of 214 time intervals. The data reduction was performed 
using the standard software (HE ADAS software, v6.1, CALDB 
version Jul07) and following the procedures reported in the in- 
strument user guide£]. The spectral analysis was performed us- 
ing XSPEC (vll.3). The 214 spectra were binned in order to 
ensure a minimum of 20 counts per energy bin, ignoring chan- 
nels below 0.3 keV and above 10 keV. For the Galactic hydro- 
gen column density A^h.mw we assumed the value reported in 
Dickey & Lockman (1990)| along the GRB direction. We mul- 
tiplied to the spectral models an extra neutral absorber at the 
source redshift letting the column density A^h„- free to vary. SPL 
provided very good fits in most cases and we found that the pho- 
ton index F ranges from 0.5 to 3.9. In particular there are 13 
spectra (6% of the total), belonging to 5 different bursts with 
F >3 . These are: one from GRB 060502A, GRB 060614 and 
GRB 060904B, three from GRB 061 1 lOA and seven from GRB 
060729. They are all collected in Windowed Timing (WT) mode 
(see Hill et al. 2004 for a description of the different operational 
modes of the XRT) and belong to the prompt-afterglow transi- 
tion phase, observed less than 500 seconds after the event trig- 
gers. In Fig. [U the 13 very soft spectra are shown together with 
their best SPL fit. The mean xl value of this 13 spectra sub- 
sample is 1.20+0.07, significantly worse than the average of the 
entire sample that is 1.01+0.01 (chance probability <10"'*). In 
particular there are some evident departures from SPL model at 
high energies in four spectra from four bursts. 



3. Excess statistical significance 

Having in mind the rules of thumb given by Protassov et al. 
(2002), to determine if there are statistically significant depar- 
tures from SPL, we followed step by step the method described 
by Rutledge & Sako (2003) which was used to calculate the sig- 
nificance of many X-ray features by |Sako et al. (2005)| First, 
we considered the redistribution matrix (RMF=RMF(PI,E)) and 
we fitted each column of the matrix (2400 in total; se^B) by a 
Gaussian function. As done by |Rutledge & Sako (2003)] we built 




' http://heasarc.nasa.gOv/docs/swift/analysis/#documentation 
^ http://heasarc.gsfc.nasa.gov/docs/heasarc/caldb/swift/ 



3.5 4.0 4.5 

Power law spectral index 

Fig. 2. Black point and black contour are the best fit and 1,2,3 cr 
confidence contour of the fit on one observed spectrum of GRB 
060729. Grey dots are the best fit results on a sample of 20,000 
simulated spectra. 

a new RMF with the 2400 columns replaced by the Gaussian fit 
of the original RMF. Because our features are broader than the 
instrumental spectral resolution (which is 0. 1 1 keV at 4 keV) we 
also built four artificial RMFs with Gaussian functions 3,5,10,16 
times wider than the best fit of the original RMF value. In prac- 
tice, we built a set of 4 different RMF worsening the spectral 
resolution to look for the scale which maximizes the signal-to- 
noise ratio. 

For each of the 13 spectra being tested, we convolved the 
PI count spectrum with the 5(1 nominal, 4 smoothed) RMFs. 
Then, for each of the 13 observed spectra, we created 100,000 
Montecarlo (MC) realizations of the raw pulse-invariant (PI) 
spectra based on the respective best-fit SPL models using a fixed 
number of source plus background photons. Background events 
were randomly selected from a background spectral model, de- 
rived from fits to spectra obtained from a source-free region of 
a deep exposure. 

To check the accuracy of our MC simulations for each GRB 
we ran the XSPEC grouping and fit procedures on a sample 
of 20,000 simulated spectra. In Fig. |2] for one spectrum of 
GRB060729 we compare the results of the SPL model fit of the 
observed data, with the results of the fit on simulated data. 

As we did for the observed data, we convolved all the 
100,000 simulated PI spectra with the 5 matrices correspond- 
ing to the 5 filter scales. For each spectrum and each scale we 
counted the number of MC realizations exceeding the feature in 
the data and recording the energy and the scale for which this 
number is minimum. There are four spectra, in the sample of 
13, for which we found excesses in the data that can be repro- 
duced as statistical fluctuations of the SPL in less than 10 tri- 
als out of 100,000. This correspond to a single trial significance 
higher than 99.99% . These are the last WT spectra, coinciding 
with the last phase of the X-ray light curve steep decay, of GRB 
060502A, 060729, 060904B and 061110A (Tableland Fig.[3] 
nil. In the rest of the paper we will limit our analysis to these four 
spectra. 

In order to improve the accuracy in the calculation of the ex- 
cess statistical significance, for the sample of four spectra we en- 
larged the simulated sample to 1,000,000. Results are reported in 
Table[T]and illustrated in Fig.|5] We note that in the case of GRB 
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Fig. 3. The spectra (in energy units) of the four GRB afterglows for which we detect significant departures from SPL model which 
is plotted in grey. 



Table 1. (I) GRB: Name of the GRB; (II) Rest-frame time 
interval; (III) Redshift; (IV) energy of the excess maximum 
(observer frame); (V) energy scale in times of instrumental 
resolution (VI) Single trial significan ce; (VII) Multi trial sig- 
nificance. Ref: '|La Parola (2006)1 jCucchiara et al. (2006) 
Grupe et al. (2007)1 ~' 



- Grupeetal. ( 2006); '^ Fugazza et al. (2006) 
iThoene et al (2006b), 



IThoene et al (2006a) 



Fox et al. (2006) 
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energy resolution scales. For each energy scale we considered a 
different number of energy resolution elements: 40,14,8,4,3 re- 
spectively for the 1,3,5,10,16 scales. The results are reported in 
Table [T] Because we searched for features on a homogeneous 
sample of 13 spectra, resulting in 4 successes, we can also es- 
timate the joint statistical significance. In particular, we set a 
lower limit to the total joint probability using the binomial dis- 
tribution and assuming that all the four spectra have the same 
multi-trial significance equal to the lowest (P=99.04%). Then 
we searched for the probability to have a rate of 4 successes out 
of 13 tests with mean probability P, resulting in a joint signif- 
icance of 99.9994% (4.3cr) which ensures that the features we 
detected cannot be explained as statistical fluctuations beyond 
any reasonable doubt. 



060729, with 1,000,000 MC tests we obtained a significance 
value slightly lower than the one we obtained with 100,000 
tests and sHghtly lower than the 99.99% threshold (99.9875%). 
Because the two results are perfectly consistent (within Icr er- 
rors) we kept GRB 060729 in our sample. To calculate the cor- 
responding multi-trial significance, we took into account that we 
searched for excesses on a sample of 13 spectra on five different 



3.1. Instrumental effects 

We investigated the possibility that these features are produced 
by some instrumental effects such as pile-up. This cannot be the 
case, because the mean count rate in all the four cases is less than 
33 counts s"' and it is always below 70 counts s"', a factor 3 be- 
low the pile-up threshold in WT (~ 200 counts s ' , see Campana 
et al. 2008). We can also exclude that these features are produced 
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Fig. 4. Upper panels: The luminosity (0.3-10 keV rest-frame energy band) curves of the four GRB afterglows for which we 
detect significant departures from SPL model. With the continuous lines we plot the best fit using the model described by 
Willingale et al. (2007), split in prompt and afterglow contributions (dashed lines). With the filled circles we plot the luminosity 
of the spectral excesses we detect. The upper limits to the detection of the extra components (see Section |4|i in the diff'erent time 
slices are also reported as top-down arrows (for the clarity of the picture we report only the value calculated assuming a DPL model, 
see SectionHJl. With the dotted line we plot the fit of the prompt emission re-normalized to the luminosity of the extra component. 
Lower panels: The spectral indexes F given by SPL fit, measured in the diff'erent time slices. 



by some anomalous hot pixels from the comparison with the ex- 
pected point spread function ( IMoretti et al. 2005T l. 

We can also exclude that the excesses are due to uncertainties 
in the instrument response calibration. All the excesses are found 
at energies > 3.5 keV (Table[T]l in the central part of the CCD (< 
70 pixels, equivalent to ~ 3 arcmin). In this position and energy 
range the systematics in the effective area, quantum efliciency 
and energy redistribution calibrations are less than 5%. (se^and 
Godet et al. 2007b). In fact all the major instrumental edges are 
below 3.5 keV (Au for the mirror, Al for the filter, C, N, O and 
Si for the CCD) and there is no evidence of position dependence 
at energies higher than 0.5 keV. 

The excesses cannot result from the background subtraction 
procedure. In fact, the expected background events, in the extrac- 
tion region we considered, in 100 second exposure, is L2+0.2; 
among them 0.6 are expected with energy higher than 2 keV 
(Moretti et al. 2007). Moreover, we checked our data against un- 
covered anomalies performing our analysis with different back- 
ground extraction regions and we found perfectly consistent re- 
sults. 

In the case of GRB060904B, a galaxy cluster with a core 
radius of 12 arc sec is present in foreground, at 2.3 arcminutes 

' http://heasarc.gsfc.nasa.gov/docs/heasarc/caldb/swift/docs/xrt 



of projected distance from the GRB. By modeling the cluster 
surface brightness with a King profile in the extraction region of 
the afterglow we expect to find <1 photon (99.99% confidence) 
from the cluster in the 80 s exposure. 

4. Modeling the data with an extra component 

As discussed in the previous section we found a sub-sample 
of four spectra, belonging to four different GRB which present 
highly significant deviations from the SPL models. In fact, if 
fitted with SPL, they all four present a clear excess at high ener- 
gies. 

As already said, most of the departures from SPL models 
found in time resolved spectral analysis of Swift observations of 
early afterglow have been explained so far by the curvature of 
the spectrum and the presence of the vFy peak (Epeait) within the 
XRT band (Falcone et al. 2006, Butler & Kocevski 2007, Goad 
et al. 2007, Mangano et al. 2007). The four spectra we are con- 
sidering can be modeled neither by a Band model (Band et al. 
1993) nor by a power law with an exponential cutoff. In the time 
intervals we are considering, the BAT signal is almost null. This 
prevented us from studying the spectrum in the combined energy 
band. However, any attempt at fitting the XRT spectra with a cut- 
off power law or with a Band model could not constrain the Epeaii 
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Fig. 5. Data are compared with the single trial confidence contours as resulted from MC simulations (10^ tests for each GRB; see 
TabllJ. In the plot, data and simulations contours are both smoothed to the energy resolution scale which maximizes the spectral 
feature signal. These are 7, 10, 3 and 5 times the instrumental spectral resolution, approximately corresponding to 0.7, 1.0, 0.3, 0.5 
keV respectively for the four afterglows (note that energy resolution varies within the XRT energy band so these are only rough 
estimates). The circles indicate where the most significant excess is detected. 



Table 2. (I)GRB: Name of the GRB; (II)Model: fitting model (see text); (III)A^h,z: neutral absorber column density at the source 
redshift; (IV)ri, (r2): photon index of the first (second, present if model is DPL) power law; second photon index is fixed to the 
value of the late afterglow; (V)E: rest frame blackbody temperature if Model is BBl, or BB2 or mean value of the Gaussian line if 
model is GAU; (VI)Ri: blackbody radius; (VlTjcTg: cr of the Gaussian line; (VIII)L: unabsorbed luminosity of the extra component 
in the rest-frame [0.3-10] keV energy band; (IX)^^(dof): reduced (degrees of freedom) in two different energy bands. 
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Fig. 6. The anomalous spectra in the GRB 060502A and GRB060729 afterglows. For the sake of clarity, we use two plots for each 
spectrum: we plot DPL (dotted grey line) and GAU (dashed line) on the left plots and BE 1 (dashed line) and BB2 (dotted grey line) 
on the right plots. In each plot we show the contribution of the single components and their sum, with the same line style and color 
SPL best models are plotted with the continuous black step-line. In the lower panels of each plot we show the ratio between data 
and models with the same color code. For the sake of clarity we plot errors only for SPL models 



within the XRT band and therefore did not improve the SPL fit. 
We found that our analysis is consistent with previous papers. 
In particular Butler & Kocevski (2007) performed time resolved 
spectroscopic analysis of BAT and XRT simultaneous observa- 
tions of a large sample of Swift GRBs. For GRB 060729 and 
GRB 060904B, in particular, they found that, during the early 
phases of the afterglow, the energy peak (Epeak) of the prompt 
emission spectrum transits in the X-ray band. The time intervals 
where we observe excesses in the spectrum of GRB 060904B 
(Table [1]) roughly correspond to the last 5 temporal bins of their 
analysis (240-3 15 seconds in the observer frame). Although they 
do not explicitly report the best fit parameter and the values, 
it is clear from their Fig. 9 (third panel from the top on the right) 
that while the Band model is a good description of the data at 
the beginning of the XRT observation, in these five particular 
intervals. Band model gives a very poor description of the data, 
always leaving one parameter unconstrained. The same conclu- 
sion can be also drawn for GRB 060729, where our time interval 
corresponds to their last three time slices. We also note that for 
this GRB the same conclusion is also confirmed by Grupe et al. 
(2007) who found that cutoff power law model gives xl values 
larger than 1 .5 in the same time interval. 

Therefore we tried to fit the data adding three different com- 
ponents to the SPL: (i) a second power law component with 
the slope frozen to the value of the late afterglow spectrum let- 
ting only the normalization vary (DPL); (ii) a blackbody; (iii) a 



Gaussian line (GAU). As it will be explained later (Sect. |5.3l l, 
we found that, with the blackbody model, two equally good fits 
could be found with quite distinct parameter values. Therefore 
we considered a blackbody with temperature varying in the 0.1- 
10 keV energy band with initial guess kT=0.2 keV (BBl) and 
kT=2 keV (BB2). Q. The results are reported in Table |2] and 
shown in (Fig.|6l|7]i. 

5. Discussion 

5.1. Energetics and time variability of tine spectrai features 

The DPL, BB2 provided an extra component to SPL which com- 
pensates the high energy residuals. In BBl models, the black- 
body component represents a significant fraction of the softer 
part of the spectrum, while the excesses at high energies are 
accounted for by the power law component. The (unabsorbed) 
luminosities of the additional components calculated in the [0.3- 
10] keV rest-frame energy band are typically 1-10% of the total, 
in DPL, BB2 and GAU models, while it is 10-20% of the total 
in BB2 models (Table|2]and (Fig.|6]|7]l). As illustrated in Fig.|4] 
for all the four afterglows we could perform time resolved spec- 
troscopy. The time intervals from which we extracted the anoma- 



* Note that with DPL, BB 1, BB2 and GAU we intend the model SPL 
plus the extra component. 
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lous spectra correspond to the last WT spectrum, right before the 
canonical steep-shallow light-curve break (Nousek et al. 20061 ). 
With the same criteria adopted for the original 13 very soft spec- 
tra we did not find any other significant deviation from SPL 
model in any of the time slices considered. In order to study 
the time variability of the spectral features, for each of the four 
GRBs, we estimated the detection upper limits. To do this, for 
each slice, we added an extra component (second power law, 
blackbody and Gaussian) with the parameter value set to the best 
DPL, BBl, BB2 and GAU values, letting the normalization free 
to vary. We set the upper limits for the detection when these ad- 
ditional components produce a factor of 3 worsening in terms 
of null hypothesis probability of the xl of the SPL fit. We note 
that this should be considered as a rough estimate of the up- 
per limits; a rigorous calculation of all the upper limits would 
have required an unrealizable number of simulations. However 
we verified that, at least in one case, the upper limits roughly 
calculated differs by less than 30% from the one rigorously cal- 
culated. 

The upper limits vary during the afterglow depending on 
the flux and on the softness of the spectrum. In particular, as 
shown in Fig.lH the detection of the spectral features, in the last 
spectrum before the plateau phase, coincides with a drop of the 
detection threshold. This is due to the simultaneous flux decay 
and spectral softening which allow the detection of spectral fea- 
tures in the higher part of the energy band. We note that GRB 
060904B did not show this component at early times, although 
the sensitivity was also good then. We also note that in at least 
three cases the extra components present in the last part of the 
light curve steep decay disappear in the shallow phase, although 



the upper limits in this time slices are very low (in the case of 
GRB 060904B, the observation stopped during the steep decay 
of a giant X-ray flares). Evidently, whatever its nature, the emis- 
sion mechanism responsible for these spectral features varies on 
a time scale similar to the prompt emission. 

5.2. Double power law (DPL) 

In three cases DPL model provides the best improvement in the 
fit, taking into account that we add only one extra parameter to 
the SPL model. This model would provide a natural explana- 
tion to the excesses we observe. In fact, because the afterglow 
spectrum is significantly harder than the prompt tails, when the 
prompt flux decreases and softens, the afterglow emission be- 
comes visible at higher energies. This would easily explain the 
fact that the excesses are detected just before the steep-shallow 
light-curve break and disappear in the following time slice. On 
the other hand, as shown in Fig. ID the excess luminosity is much 
larger than the expected contribution of the afterglow forward 
shock component alone (e.g. lSari 1997 ', Willingale et al. 2007). 
But it might be explained by the radiation produced by an ex- 
treme reverse shock in the X-ray band (see Zhang et al. 2006, 
Kobayashi & Zhang 2007). 

5.3. Power law plus blackbody (BB1, BB2) 

Good fits are also provided by adding to the SPL a blackbody 
component (2 extra parameters). In all the four afterglows, with 
this model, the fit results depend on the blackbody tempera- 
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Fig. 8. Confidence contours (68% and 90%) of the single power 
law+blackbody model for the four different GRBs. 



ture initial guess. As shown in Fig. [8] has two (and only 
two) different but equally significant minima. BBl, with initial 
guess kT=0.2 keV, gives good fits with (redshift-corrected) tem- 
peratures in the range 0.20-0.33 keV, radii in the range (4.9- 
7.3)xl0'^ cm and power law indexes in the range 2.2-2.9. BB2, 
with initial guess kT=2 keV gives good fits with red-shift cor- 
rected temperatures in the range 1.3-3.2 keV, radii in the range 
(1.4-5.6)xl0"^ cm and power law indexes in the range 3.8- 
4.5. For GRB 060502 A and GRB 06111 OA, BBl fits provided 
slightly better as calculated on the whole energy band, while 
for GRB 060729 and GRB 060904B the fits yield equivalent re- 
sults. Interestingly, in each of the four cases, the;^^^ matrix pro- 
jection on the photon index - temperature parameter plane has 
always two well defined minima which are split in two different 
regions of the plane (Fig.[8]l. This explains why fit results depend 
on the initial guess value of the blackbody temperature. 

In the case of GRB 060729 our BB 1 model is consistent with 
the results of Grupe et al. (2007) and Godet (2007a). They in- 
terpreted the thermal emission as due to the photospheric emis- 
sion from X-ray flares. In fact, assuming that X-ray flares, in 
the early afterglow, are produced by the same mechanism of the 
prompt phase, a thermal component in the early afterglow spec- 
tra is expected in a similar way to the prompt emission (Ryde et 
al 2006). This interpretation can be easily extended to the other 
three GRBs that have very similar characteristics. The prompt 
thermal emission described by BB2 model can also be explained 
as a shock break-out. The shock-heated plasma would be at 
temperatures of (1.4-3.7) x 10^ K, with a luminosity of (1.0- 
7.0)xl0'*^ erg s"' corresponding to a radius of (1.4-5. 6)xl0'" 
cm. Assuming the duration of the time slice in which we ex- 
tracted the spectrum as the duration of the emission, we obtain 
a total energy of (1.0-1.8)xlO'*' erg. The energy, variability, lu- 
minosity and temperature we observe in the detected excesses 
are consistent with the characteristics of the transient event from 
shock breakout in Type Ibc supernovae, produced by the core- 
collapse of WR stars surrounded by dense winds dLi 20071 ). 



5.4. Power law plus a Gaussian line (GAU) 

Adding a Gaussian line to the SPL in three cases does not im- 
prove the fit with respect to BB2, although it uses one extra pa- 
rameter. In fact, for GRB060502A, 060729 and 061110A the 
best fit is given by low energy and very broad lines (1-3 keV 
with cr = 2-8 keV), with best fit values poorly constrained. The 
case of GRB060904B is different and much more intriguing: 
here the Gaussian fit provides very well constrained values for 
the line. In the rest frame the mean value is 7.85^^25 keV, the 
widfli 0.50:!:U-5^ keV and its luminosity is (1.10+0.18)x lO'^^erg 
s"'. Interestingly, the Gaussian component can be explained as a 
line emission of highly ionized Nickel (7.81 keV). We refer the 
reader to an accompanying paper (Margutti et al. 2007) for a de- 
tailed discussion of the theoretical implications of the possible 
detection of Nickel emission at ~ 200 sec after the onset of the 
GRB. 



5.5. A Nickel line in the GRB 060904B afterglow spectrum ? 

As we saw in the previous sections, the DPL, BBl and BB2 
models provided significant improvement in the fits with respect 
to SPL model in all four cases. In the GRB 060904B afterglow 
spectrum these models left some residuals in the high energy 
part of the spectrum which can be fit well only by GAU model 
(Fig-0- Since we cannot lean only onx^ statistics to evaluate the 
goodness of the DPL, BBl, BB2 fits, we tested the probability 
that these residuals are statistical fluctuations of the DPL, BBl, 
BB2 models. 

To this aim, we adopted the same procedure we previously 
used to test the hypothesis that the residuals were fluctuations of 
a SPL model (Sect.[3]l. We only replaced the SPL model with the 
DPL, BBl, BB2 as input model for the MC simulations. When 
we tested the GRB 060904B afterglow residual significance as a 
fluctuation of a SPL model, we found 4.2 cr (i.e. 99.9993%) as 
a single trial, corresponding to 3.2 cr (i.e. 99.9461%) as multi- 
trial (see Table [U. If we test, instead, the possibility that this is 
a statistical fluctuation of a more complex spectral model (DPL, 
BBl, BB2) the single (multi-) trial statistical significance of this 
detection is 2.7 (2.2) cr, with very small differences among the 
three models. This means that the deviation from SPL model that 
we observe in GRB 060904B can be described as a Gaussian 
deviation from a SPL model at 3.2cr or from a two component 
model at 2.2 cr. 

6. Conclusion 

Our most solid result is that we found a small homogeneous and 
fairly defined sample of afterglow spectra (the soft sample) for 
which deviations from the SPL spectral model are highly proba- 
ble. We started from an homogeneous sample of bright GRB af- 
terglows with known redshift and we studied their spectral evo- 
lution. We split the data in different time slices and we focused 
on the softest spectra. In this sub-sample at least 4 cases out 
of total of 13 present highly significant deviations from the SPL 
spectral model during the prompt-afterglow transition phase. We 
could firmly exclude that these excesses can be explained as sta- 
tistical fluctuations of a SPL spectrum or some instrumental ef- 
fects. We also excluded that data can be fitted by a Band or a 
cutoff power law model. We fitted these spectra adding one of 
three different trial components to the SPL model. We did not 
try to discriminate among these different models on a purely sta- 
tistical basis, and we discussed them using the component time 
variability and energetics. 
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In a very recent paper Yonetoku et al. (2007) show that 
the very soft spectrum of the early afterglow of GRB 060904A 
presents a feature which is very similar to what we described in 
the present work. They selected and stacked the data in the time 
intervals of the GRB 060904A early afterglow, where the spec- 
tral photon index is larger than 4.0 (we note that this GRB is not 
included in our sample because its redshift is not known). In a 
very similar way to our results, they found that, in this spectrum, 
the data show a clear hardening break around 2 keV. This feature 
leaves a significant excess with respect to the double power law 
model above 2 keV. They conclude that this spectrum consists of 
two emission components, the second being consistent with the 
spectrum of the late afterglow (our DPL model). These two (in- 
dependent) studies represent a direct piece of evidence that the 
emission observed in the early phases of the afterglow is com- 
posed by more than one component. Differently from Yonetoku 
et al. (2007), in our sample, we showed that, if the second com- 
ponent were the emerging afterglow emission, at early time this 
should be much more luminous than the expectation from the 
classical afterglow model. We showed that spectral studies of 
the prompt-afterglow transition phase can be the starting point 
to seperate different emitting components and could provide use- 
ful information in order to better understand the afterglow light 
curve complexity. 
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